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ABSTRACT 

Formation of monodispersed colloidal particles is a complex process: 
nuclei, produced rapidly in a supersaturated solution, grow to nanosize pri- 
mary particles, which then aggregate (coagulate) to form much larger final 
colloids. This paper reviews a kinetic model that explains the formation of 
dispersions of narrow size distribution in such processes. Numerical simu- 
lations of the kinetic equations, with experimental model parameter values, 
are reported. The model was tested for a system involving formation of 
uniform spherical gold particles. The calculated average size, the width 
of the particle size distribution, and the time scale of the process, agreed 
semiquantitatively with the experimental values. 

1. INTRODUCTION 

We review the recently developed theoretical model [1] of growth of 
monodispersed colloids by precipitation from homogeneous solutions [2-8]. 
Systematic experimental studies of the mechanisms of synthesis and prop- 
erties of such colloids have been initiated about quarter a century ago. A 
large number of uniform dispersions of particles of various chemical com- 
position and shape, ranging in size from nanometers to few microns, have 
been described. Theoretical description of the formation of uniform colloids 
relies on experimental identification of solute species that are involved in 
the various stages of the particle formation and on observation of growth 
stages leading from the initial nucleation to final particles. Early theoret- 
ical modeling was based on the mechanism suggested by LaMer: a short 
nucleation burst, followed by diffusional growth of the resulting nuclei to 
form identical fine particles [8,9]. 

However, there has been mounting experimental evidence that the 
burst-nucleation/diffusional growth mechanism alone is inadequate. Specif- 
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ically, it has been found that many spherical particles precipitated from so- 
lution showed polycrystalline X-ray characteristics, such as ZnS [10], CdS 
[11], Fe203 [12], etc. These particles are not single crystals. Rather, it has 
been confirmed by several techniques (small angle light scattering, electron 
microscopy, X-ray diffraction) that most monodispersed colloids consist of 
small crystalline subunits [10-20]. Furthermore, it has been observed [1,7,16] 
that the crystalline subunits in the final particles were of the same size as 
the diameter of the precursor subunits (singlets) formed in solution, thus 
suggesting an aggregation-of-subunits mechanism. 

The substructure have been identified also in particles of shapes differ- 
ent from spherical, and it has been recognized that different morphologies 
of the final colloids must be related to the nature of the precursor singlets 
[12,18-20]. This experimental evidence has led to two major theoretical chal- 
lenges. First, the morphology and shape selection of particles formed by 
interplay of nucleation and aggregation processes must be explained. Sec- 
ond, the size-selection mechanism, i.e., the kinetics of generation of narrow 
particle size distribution, must be identified. Several theoretical approaches 
utilizing thermodynamic and dynamical growth mechanisms [1,4,5,18,21-32] 
have been described in the literature. Models of aggregation of subunits can 
be developed that yield a peaked and even sharpening with time particle 
size distribution [21-34]. However, none of the earlier attempts could fit 
quantitatively a broad range of experimental findings. Here we review the 
first results of a new approach [1] that explains the size selection, by cou- 
pling the dynamical equations for the processes of the particle growth by 
aggregation of subunits and of formation (nucleation) of these subunits. 

Thus, the main new finding of our work [1], crucial to obtaining narrow 
size distribution, has been that the growth of the final, secondary particles 
by aggregation of singlets, must be coupled, via the rate of formation of 



- 3 - 



these primary particles (singlets), to the time-dependence of the process 
of their nucleation and growth. We take the simplest possible models of 
both processes (primary and secondary particle formation). This choice 
simplifies numerical simulations and thus allows to scan a wider range of 
parameters. It avoids introduction of unknown microscopic parameters; as 
a result we only fit one parameter, the effective surface tension, and even 
that parameter turns out to be close to the experimental bulk value. This 
review is organized as follows. In Section 2, a kinetic model of secondary 
particle formation by singlet-capture dominated growth is considered. In 
Section 3, primary particle (singlet) formation by burst nucleation is in- 
corporated to complete the model. Finally, Section 4 reports results of 
numerical simulations and comparison with an experiment. 

2. KINETIC MODEL OF COLLOID PARTICLE GROWTH 

In colloid particle synthesis of the type considered here, in the ini- 
tial induction step solutes are formed to yield a supersaturated solution, 
leading to nucleation. The nuclei might then further grow by a diffusive 
mechanism. The resulting primary particles (singlets) in turn aggregate to 
form secondary particles. The latter process, to be modeled in this sec- 
tion, is facilitated by the appropriate chemical conditions in the system: 
the ionic strength and/or the pH must assume values such that the surface 
potential approaches the isoelectric point, resulting in reduction of electro- 
static barriers, thus promoting particle aggregation. Formation of the final 
(secondary) particles, which can be of narrow size distribution, is clearly a 
diffusion-controlled process [1-7]. 

An important experimentally observed property is that the secondary 
particles are sufficiently sparsely positioned in solution to consider their 
evolution as largely independent. Furthermore, we will assume that the 
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particles are spherical, with density close to that of the bulk material, which 
is experimentally a very common case, and the focus will be on their size 
distribution due to growth by aggregation, essentially, by consumption of 
diffusing singlets (monomers, primary particles). Several other simplifying 
assumptions will be made in order to zero-in on the essential ingredients 
of the theoretical modeling and make numerical simulations tractable. For 
instance, we do not address the processes by which the singlets forming 
secondary particles undergo restructuring and rearrangements resulting in 
compact structures and leading to shape selection. Some of the approxi- 
mations made will be discussed at the end of this section. Thus, we focus 
on the two main stages in the process: the production of the primary and 
secondary particles. The latter will be treated presently, while the forma- 
tion of nuclei and their possible growth (aging) to primary particles will be 
considered in the next section. 

Kinetic rate-equation models of aggregation typically utilize a master 
equation for the distribution of growing particles by their size. Here the size 
will be defined by how many primary particles (singlets) were aggregated 
into each secondary particle, denoted by s = 1, 2, . . .. The growing particles 
can adsorb or emit singlets and multiplets. The master equation is then 
quite standard to set up, as in, e.g., [33]. In the present case, the process 
is experimentally documented to be highly irreversible, so that detachment 
can be disregarded. By considering the singlet number s as the only param- 
eter characterizing the cluster, we use the approximation that assumes that 
the internal relaxation/restructuring processes within the cluster (growing 
secondary particle) are fast, and that they result in an approximately spher- 
ical object, of density similar to that of the bulk material, so that the voids 
have been eliminated and the solvent pushed out in the internal restructur- 
ing. We will comment on this point again later. 
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Furthermore, it is assumed that the singlets have diffusion constants 
larger than aggregates so that their capture dominates the growth process; 
this approximation will be also discussed again, shortly. The master equa- 
tion then takes the form 

d ^ = w s _ 1 N s _ 1 -w s N s (s>l). (1) 

Here N s (t) is the time-dependent number density (per unit volume) of the 
secondary particles consisting of s primary particles. The attachment rate 
w s will be taken from the Smoluchowski expression [35]: 

w s (t) = AnRsDN^t) , (2) 

where R s is the radius of the s-size particle, given by 

R s = 1.2 rs 1/3 . (3) 

The parameters r, the primary particle radius, and D, their diffusion con- 
stant, are experimentally available. The constant 1.2 in equation (3) was 
calculated as 

(0.58)~ 1/3 ~ 1.2, (4) 

where 0.58 is the typical filling factor of the random loose packing of spheres 
[36] . The fact that all the primary particles were assumed to have the same 
radius r will be further commented on later. Note that the rate in equa- 
tion (2) involves the product R S D oc rD, which, according to the Einstein 
formula for the diffusion constant, is not sensitive to the distribution of 
values of the radii r. 

We note that this approach already involves various assumptions. For 
example, to derive equation (2), the attachment process is assumed to be 
fully irreversible. Furthermore, the above equations are only valid for large 
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s. Indeed, in reality clusters of all sizes s and S can coagulate to form larger 
clusters of size s + S. These are produced with the rate ^tt{R s + Rs)(D s + 
Ds)N s Ns within the Smoluchowski approach, where for s = S there is an 
additional factor of 1/2 owing to double-counting. Our ignoring all but 
the singlet capture processes, S = 1, is supported by the experimental 
observations but constitutes an approximation. Furthermore, the form (2) 
is only correct for s ^> S = 1. Thus, we used the above relations as shown 
simply to avoid introducing additional parameters. For instance, for s = 1, 
this ammounted to ignoring the "double counting" factor of 1/2, and the 
added factor of 4 due to the radii and diffusion constants of the two singlets 
adding up. 

Ordinarily, the evolution of the population of singlets, which is not 
covered by equation (1), is obtained from the conservation of matter: 



which assumes that initially, at time t = 0, there are only singlets. When 
combined with the assumption of the singlet-capture dominance, the par- 
ticle distribution is confined to small sizes, as has been confirmed by nu- 
merical simulations and other studies [1,37]. Indeed, most singlets will 
simply combine into dimers, fewer trimers, etc., and then the growth stops. 
The conventional approach has been to consider more general models, with 
discrete or continuous population balance, restoring multiplet aggregation, 
etc., which adds terms in and modifies the master equation (1); see, e.g., 
[21-34]. Typically such models have yielded wide size distributions [33,38- 



We have developed [1] a new approach based on the observation that 
the supply of singlets is in itself a dynamical process. Numerical simula- 
tions indicate that if the concentration of the singlets were constant, i.e., if 




(5) 



40]. 



-7- 



they were continuously generated to compensate for their depletion due to 
aggregation, then the resulting particle size distribution would be wide and 



cut-off s-value. However, if the supply of singlets is controlled in such a 
way that their concentration is a decaying function of time, then size selec- 
tion can be obtained for the secondary particles with some time-dependence 
protocols. For example, when the rate pit) at which the primary particles 
are formed (per unit volume) was chosen to be a decaying power-law, then 
numerical calculations [1] yielded a single-hump size distribution for the 
secondary particles. Note that the equation for N\{t) must be modified by 
replacing relation (5) by 



with the initial values N s (0) = for all s = 1, 2, 3, . . .. 

Equations of this type, with singlet-capture dominance of the dynamics 
and several ad hoc singlet-input rate functions, have been described in the 
literature [37] . The emphasis [37] has been on cases which are exactly solv- 
able, e.g., one-dimensional versions, and those where the size distribution 
shows self-similar behavior. Our approach is quite different in that we actu- 
ally model the primary-particle input rate pit) by using the burst-nucleation 
approach; see the next section. 

Let us now further comment on the approximations involved in using 
the simplest kinetic equations for the aggregation process. These include, 
for instance, ignoring multiplet mobility and multiplet-miltiplet collisions, 
as well as the effect of mobility of aggregates (multiplets) on the diffusion 
constants used in the rate expressions, etc., especially in the beginning of the 
process when most aggregates are small. There are established methods in 



peaked at small sizes, such that N s ~ N\s x / 3 for sizes up to some growing 
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the literature that avoid some of such difficulties [21-34] and, in fact, several 
of the models for various reaction rates lead to particle-size distributions 
peaked and even sharpening with time [22,34,42,43]. However, the main 
point of our work, which we believe is new and crucial to obtaining narrow 
size distributions, has been that the growth of the secondary particles must 
be coupled, via the rate of generation of singlets (primary particles), to the 
time-dependence of the process of formation of the latter; see the following 
sections. 

Thus, we intentionally took the simplest possible models of both pro- 
cesses, the primary and secondary particle formation. This choice has the 
following advantages: it simplifies numerical simulations and thus allows to 
scan a wider range of model parameters; it avoids introduction of unknown 
microscopic parameters. As a result, for instance, our description of the 
aggregation process in this section has no adjustable parameters; they are 
all experimentally available. Similarly in the following sections, for primary 
particle formation, we only utilize one adjustable parameter, the effective 
surface tension, and even that parameter turns out to be close to the ex- 
perimental bulk value. 

However, we recognize that more sophisticated modeling can improve 
consistency with experiment, perhaps at the expense of additional assump- 
tions and parameters, and we intend to explore this avenue of investigation 
in future studies. In fact, we carried out preliminary numerical simulations 
allowing for the dimer diffusion and attachment to larger aggregates and 
showing trend of improved consistency with the experiment. We further- 
more note that the approximation of restricting the aggregation process 
to only the smallest particles sticking to the larger particles has been al- 
ready used in the literature, e.g., [22,27]. Considerations of colloid stability 
have been typically utilized to justify such approximations, and we note 
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that detailed arguments of this sort would require additional microscopic 
parameters in the model. 

We also comment that generally in diffusion limited growth the aggre- 
gates are expected to be fractal [38-41]. In this work we avoid the issue of 
shape selection and internal structure of the secondary particles; see a re- 
view [21]. There are likely some internal rearrangements and pushing out of 
the solvent going on during the secondary particle growth for typical exper- 
imental conditions, because the particles are clearly spherical throughout 
the growth process. In some experiments [14] the surface of the secondary 
particles was initially "hairy" and it got smoother at later times. We have 
assumed that the internal rearrangement processes are fast enough so that 
the shape and morphology of the aggregates are, respectively, spherical and 
compact. One could propose that the effects of internal rearrangements in 
our model make the coefficient in equation (3) another fit parameter rather 
than a fixed number. Distribution of the primary particle radii could also 
affect the porosity properties and thus modify this coefficient; we have not 
explored this matter. 

3. NUCLEATION AND GROWTH OF THE PRIMARY 
PARTICLES 

In this section we evaluate the primary-particle production rate, p(t), 
assuming a fast nucleation process [8,9,34,44]. Thermodynamic models of 
"burst" nucleation of uniform dispersions have been originally formulated in 
[8,9]. Modeling the rate of formation of primary particles (singlets) requires 
setting up a master equation, where the rate of growth is determined by 
the Boltzmann factor with the thermodynamic free energy difference AG, 
multiplied by —l/(kT), in the exponent. This approach in turn requires 
modeling of the free energy of the growing embryos (sub-critical nuclei); in 
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the simplest approach one can use the generic volume-plus-surface energy 
expressions. 

Let us refer to the species (atoms, ions, molecules) which serve as 
monomers for the primary-particle nucleation as solutes. For a given con- 
centration c(t) of solutes, larger than their equilibrium saturation concen- 
tration Co and approaching Co for large times t, the rate of formation of 
critical nuclei can be written as [34,44] 



which is based on the diffusional capture of solutes, whose effective radius 
is denoted by a, diffusion constant by T>, and n is the number of solutes 
in an embryo. The subscript ens refers to values calculated at the critical 
nucleus size. Note that c(t) = n\{t). 

The expression (7) involves the following assumptions. For embryos 
of size n < n cns , the solutes can be captured and emitted fast enough so 
that the size distribution is given by the equilibrium form. Thus, the factor 
ce - AG cns/kT ^ n e q Ua tion (7) follows from the expectation that embryo sizes 
up to n cns are thermodynamically distributed, according the Boltzmann 
form. For sizes larger than n cns , the dynamics is assumed to be fully ir- 
reversible, corresponding to an unbound growth by the capture of solutes. 

1 /3 

The factor <iiran c !n S Vc in equation (7) is thus the appropriate version of 
the Smoluchowski growth rate similar to that in equations (2) and (3). 
The filling-fraction correction factor was absorbed in the definition of the 
effective solute radius a to simplify the notation; it will be specified later. 

For the free energy of the n-solute embryo, the following expression 
will be used: 



p(t) = 47ranJ/ 3 s Pc 2 e- AG -/ fcr 



(7) 



G = -nkT In (c/c ) + 47ra 2 n 2/3 



(7 



(8) 
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which involves the bulk term, proportional to n, and the surface term. The 
standard form of the bulk term was derived as follows. It is assumed that 
the entropic part of the free-energy change between the solid and solution 
phases can be calculated as the entropy in the supersaturated liquid suspen- 
sion of solutes of concentration c, as given by the dilute (noninteracting) 
expression of the "entropy of mixing," defined, e.g., in [45]. The surface 
term in equation (8) corresponds to the assumption that the growing em- 
bryos are spherical, of radius an 1 / 3 , and introduces their effective surface 
tension a, which is usually assumed to be comparable to the bulk surface 
tension. It will be shown later that in the present case the results are very 
sensitive to the value of a. 

Obviously, all the above expressions are only meaningful for large n. It 
has been a common practice in the literature to use them for all n, as one 
of the approximations involved in a model. In what follows, we will in fact 
ignore the difference between G and AG = G(n) — G(l), and occasionally 
treat n as a continuous nonnegative variable. Both n cns and AG cns are 
explicit functions of c(t), 



T^cns — 



87ra 2 cr 



3kTln (c/cq) 



(9) 



ag = 256AV3 do) 

° ns 27(/cr)2[ln(c/co)] 2 ' 1 } 

where the critical value n cns was calculated from dG/dn = 0. 

Next, we account for the decrease in the concentration of solutes ow- 
ing to the formation of critical nuclei. Ordinarily, for n > n cns the primary 
particles grow (age) largely by absorbing diffusing solutes, and as in the pre- 
ceding section we ignore here more complicated processes such as capture 
of small embryos, dissolution, etc. Simultaneously, the primary particles 
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are also captured by the secondary particles. In the present model, it is 
assumed for simplicity that the primary particles are captured fast enough 
by the growing secondary particles so that the effect of their aging on the 
concentration of solutes can be ignored. Furthermore, it has been gen- 
erally recognized that aging, when significant, tends to sharpen the size 
distribution [21,34]. Thus, the primary particle radius r, introduced in the 
preceding section, will be assumed to have a single, experimentally deter- 
mined value, although in reality [1] they have a finite-width, albeit not very 
wide, size distribution. This approximation was already commented on ear- 
lier; it works largely because only the product rD matters in the rates in 
equations (1) through (3). 
Thus, we write 

dc 

-j t = -n cns p{t) , (11) 

which means that the concentration of solutes is "lost" solely due to the 
irreversible formation of the critical-size nuclei. Collecting all the above 
expressions, one gets the following equations for c(t) and pit): 

dc _ 163847r 5 aV 4 £>c 2 J 2567r 3 aV 3 
dt ~ ~81(/cT) 4 [ln (c/c )] 4 ^ 1 ~ 27{kTf [In (c/c )f ( ' 



, _ 327T 2 a 3 aVc 2 J 2567r 3 aV 3 , 

m - 3/cT In (c/c ) eXP i "27(/cT)3[ln(c/c )] 2 ' ' ' } 



It should be noted that replacing the distribution of the primary particle 
sizes by a single, experimentally measured average value of r, and ignoring 
their growth (ageing) after they achieve the critical size but before they are 
captured, violates the conservation of matter. Thus, in the present model 
only the shape of the secondary particle size distribution is relevant. The 
absolute number densities N s must be rescaled to correspond to the actual 
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amount of the solid matter per unit volume. The latter data are usually 
available experimentally. 

4. RESULTS AND COMPARISON WITH AN EXPERIMENT 

In order to test the model developed in Sections 2 and 3, we used ex- 
perimental data on the growth of dispersions of submicron spherical gold 
particles [1,7], which were produced by the reduction of chloroauric acid 
(HAuCU) with ascorbic acid. The simplicity of this system and the pos- 
sibility to either measure or estimate all the necessary parameters make it 
ideally suited for testing the theoretical model. In synthesis in concentrated 
dispersions, the aggregation process is assured, resulting in the formation 
of particles consisting of a large number of subunits. The precipitation pro- 
cedure used in [1,7] has resulted in spherical gold particles of a narrow size 
distribution. After a short induction period (up to 6-8 sec) nucleation oc- 
curs, followed by immediate aggregation. The total process time varies from 
3 to 20 sec, depending on the experimental conditions selected. Scanning 
electron micrographs of the resulting particles can be seen in [1,7]. 

The field emission microscopy [1,7] revealed the presence of the sub- 
units, having an approximate size of 30 to 40 nm. The calculated value for 
the packing fraction of the subunits in the aggregated secondary particles 
was characteristic of a random loose packing, usually expected in the for- 
mation of rapidly assembled systems [36] . The size of the primary particles 
was estimated from the X-ray diffraction measurements. Calculations have 
generated values between 30 and 42 nm, which were in excellent agreement 
with the subunit size data from electron microscopy. Information on the 
size distribution of the primary particles is also available [1,7]. 

Numerical simulations required to follow the time evolution of the ki- 
netic equations turned out to be large-scale. Therefore, the testing of the 
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model was restricted to one randomly selected set of experimental parame- 
ters. However, we actually varied numerically all the parameter values and 
found that the calculation results are affected to various degree by them. 
The parameter values will be discussed roughly in the order of increasing 
sensitivity of the numerical results. For consistency all input data are in 
the MKS system of units. 

The radius of the primary particles, 

r = 2.10 • l(T 8 m, (14) 

was obtained experimentally, as described in the preceding section, which 
is within the range of 0.5 • 10 -8 m to 5 • 10 -8 m, typical for the system under 
consideration. The value in equation (14) was for the experiment [1] for 
which the initial concentration c(0) was 

c(0) = 6.0-10 25 m- 3 ; (15) 

c(0) was calculated from the concentration of the gold solution used in the 
preparation of the dispersion, and it yields the initial condition in equations 
(12)-(13). The diffusion constant D of the primary particles was obtained 
from the Einstein formula, 

D = 1.03 •10" 11 m 2 sec" 1 , (16) 

with 

/cT = 4.04 • KT 21 J . (17) 

The saturation concentration of gold in solution, Co, is not well known 
and is expected to depend somewhat on the experimental conditions. Using 
2 • 10 -12 mol dm -3 [46] yields 

c ~ 1 • 10 15 m- 3 . (18) 
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The results for the particle size distribution are not particularly sensitive 
to this parameter, because it enters under logarithm in equations (12)-(13). 
The solute diffusion constant, 

P = 1.5-10- 9 m 2 sec- 1 , (19) 

was estimated similarly to D in equation (16), using the Einstein formula, 
with the radius of the gold atom of 1.44 • 10 -10 m [47,48]. The applicability 
of this formula to single atoms may not be exact, but the particle size 
distribution is not too sensitive to this parameter value: a decrease in V 
shifts the calculated distribution to somewhat smaller aggregate sizes. 

It was established numerically that the size distribution of the sec- 
ondary particles was sensitive to the values of a, the effective atomic (so- 
lute) radius, and to the surface tension a. Note that a was defined to relate 
the number of solutes n in a growing primary particle to its radius, given 
by an 1 / 3 . It is assumed that the primary particles are largely crystalline; 
thus, the best choice of a is such that 47ra 3 /3 is the volume per atom, in- 
cluding the attributable part of the surrounding void volume, in bulk gold. 
Consequently, 

a= 1.59-l(T 10 m (20) 

was obtained by dividing the radius of the gold atom (1.44- 10 -10 m) by the 
cubic root of the volume filling fraction, 0.74, of the crystalline structure of 
gold [36]. 

The effective surface tension of nanosize gold embryos in solution, cr, 
profoundly affects the numerical results. Unfortunately, even the bulk-gold 
value, which is of order 

a~ 0.58 to 1.02 N/m, (21) 
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is not well known [49], and it may differ from that of the nanosize solids. 
Given this fact, a was chosen as the only adjustable parameter in the model. 

Experimentally, the time scale on which the secondary particle growth 
effectively terminated was about 8 to 10 sec, which does not include the 
"induction" stage. In Figures 1 though 3, the results of the numerical sim- 
ulations of the kinetic equations are presented with parameters as specified 
above, for three different values of a, which clearly demonstrate the sensi- 
tivity to the choice of this parameter. In Figure 1, the case a = 0.51 N/m 
illustrates growth that already reached saturation for times up to 10 sec. 
It should be noted that the distribution evolves quite slowly with time. 
Initially, it is heavily weighed in the small-aggregate regime. Later on, 
the large-size peak develops and eventually dominates the distribution. By 
varying a near the expected range, given in equation (21), it was found 
that, for times up to 10 sec, all a values yielded smaller average sizes than 
the experimentally measured one, 

(^)average (experimental) = L0 ± 0.1 /mi . (22) 

Seeking a that would give the largest secondary particle size resulted in 

(^fitted = 0-57 ±0.04 N/m, (23) 

which agrees well with the bulk value in equation (21). 

Figure 2 shows the size distribution for a = 0.57 N/m. The growth did 
not reach the full saturation at the relevant times, and the peak particle 
radius at t = 10 sec, of R s ~ 0.32 //m, is smaller than the experimental 
value in equation (22). The width of the distribution, of ~10%, is close to 
that established experimentally. Given the approximations involved in the 
model, only semiquantitative agreement with the experimental data should 
be expected. Since the key feature of the model is the prediction of the 
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narrow- width distribution of secondary particle sizes, the overall consistency 
with the experimental results is gratifying. 

As the value of a is increased, the large-size peak does not fully develop 
on the relevant time scales, as exemplified by the case a = 0.63 N/m shown 
in Figure 3. The reader should be reminded that, owing to the absence 
of the conservation of matter in this model, the number densities N s must 
be rescaled according to the actual amount of the solid matter per unit 
volume, if the comparison of the calculated and experimental distribution 
is attempted. 

Asymptotically, the particle-size distribution "freezes" for large times, 
i.e., the particle growth actually stops in this model as opposed to the self- 
similar growth studied, for instance, in [37]. Equation (11), with (9), can 
be integrated in closed form to yield 



where 




The left-hand side of equation (24) is the total number of primary particles 
produced by the time t. From equations (24)- (25), this number is finite as 
t — > oo, when c{t) — > Cq. The supply of the primary particles, manifested 
by the peak at small sizes for short times, in Figures 1-3, is essential for the 
large-size peak in the distribution to develop and grow, because the present 
model assumes the growth of the secondary particles to be solely by singlet 
capture. 

Let us consider the quantity r defined [1] by 

t 

T = J Nxit^dt' . (26) 
o 



If the independent variable is changed from t to r in equations (l)-(4), 
the resulting relations are linear in N s> \. One can establish that the only 
way to have a normalizable stationary size-distribution as t — > oo is to 
have r(t) approach a finite value for large t. This quantity was calculated 
numerically for the a values used in Figures 1-3. The results, shown in 
Figure 4, confirm the earlier observation that the growth process saturates 
fast for a = 0.51 N/m. For the two larger a values there is still some 
variation for the time scales of order 1 to 10 sec; see Figure 4. The function 
r(t) is useful in identifying the time scales of the growth process. 

In summary, we reviewed a new model explaining synthesis of submi- 
cron size polycristalline colloid particles with the size distribution that is 
peaked at an average value corresponding to a large number of primary 
particles in a final secondary particle. For the experimental gold-sol sys- 
tem the model has worked reasonably well: the average size, the width of 
the distribution, the time scale of the process, and even the fitted effective 
surface tension were all semiquantitatively consistent with the measured or 
expected values. The present model has involved several simplifying as- 
sumptions. Future studies will incorporate additional effects in the model 
and test it for a wider range of experimental systems. The main conclusion 
has been that multistage growth models can yield size-selection as a kinetic 
phenomenon, which has been observed in a large number of experimental 
systems. 
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Figure 1: Distribution of the secondary particles by their sizes, calculated 
for times t = 0.001, 0.01, 0.1, 1, 10 sec, using a = 0.51 N/m. 
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Figure 2: The same plot as in Figure 1, but using a = 0.57 N/m, for times 
t = 0.1, 1, 10 sec. 
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Figure 3: The same plot as in Figure 1, but using a = 0.63 N/m, for times 
t = 0.1, 1, 10 sec. 
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Figure 4: The function r{t) calculated for the values of a corresponding 
to Figures 1-3. 

-26- 



